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Abstract 

Darwin's naturalization liypothesis (DNH) proposes tliat colonization is less lil<ely wlien tlie colonizing species is related to 
members of the invaded community, because evolutionary closeness intensifies competition among species that share 
similar resources. Studies that have evaluated DNH from correlational evidence have yielded controversial results with 
respect to its occurrence and generality. In the present study we carried out a set of manipulative experiments in which we 
controlled the phylogenetic relatedness of one colonizing species [Lactuca sativa) with five assemblages of plants (the 
recipient communities), and evaluated the colonizing success using five indicators (germination, growth, flowering, survival, 
and recruitment). The evolutionary relatedness was calculated as the mean phylogenetic distance between Lactuca and the 
members of each assemblage (IVIPD) and by the mean phylogenetic distance to the nearest neighbor (MNND). The results 
showed that the colonization success of Lactuca was not affected by MPD or MNND values, findings that do not support 
DNH. These results disagree with experimental studies made with communities of microorganisms, which show an inverse 
relation between colonization success and phylogenetic distances. We suggest that these discrepancies may be due to the 
high phylogenetic distance used, since in our experiments the colonizing species {Lactuca) was a distant relative of the 
assemblage members, while in the other studies the colonizing taxa have been related at the congeneric and conspecific 
levels. We suggest that under field conditions the phylogenetic distance is a weak predictor of competition, and it has a 
limited role in determining colonization success, contrary to prediction of the DNH. More experimental studies are needed 
to establish the importance of phylogenetic distance between colonizing species and invaded community on colonization 
success. 
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introduction 

Biological invasions have attracted the attention of modem 
ecologists and biogeographers [1] because of their leading role as 
components of global change [2] . At present, organisms belongmg 
to diverse taxonomic groups are being translocated from one 
region to another with which they do not share an evolutionary 
history [3,4] . Although it is estimated that most of the organisms 
that start this dispersion do not get to become established 
successfully, sometimes they can constitute a founding colony 
and become naturalized [5,6]. One of the central challenges in the 
study of biological invasions has been to understand what factors 
determine this naturalization process [3-4], [7]. 

Various hypotheses have been proposed to explain why some 
species are capable of colonizing successfully (i.e., become 
naturalized) whUe others are not [6], [8-9]. A particularly 
intriguing and controversial role is that played by Darwin's 
naturalization hypothesis (DNH), which states that naturalization 
success depends on the phylogenetic relatedness between the 



colonizer and the members of the recipient community [10]. If this 
relatedness is close, then the colonization process will be inhibited 
as a result of the greater competitive intensity that there is — 
supposedly- between closely related species [1 1]. Conversely, if the 
phylogenetic relationship is distant, the establishment would be 
favored as a result of a lower competitive intensity. Underlying this 
relationship between phylogenetic distance and invasion success, it 
is assumed that closely related species shares similar resources and 
natural enemies [11]. 

Following the influential paper by Daehler [12], DNH has 
received renewed interest, generating controversy with respect to 
its explanatory value on the invasion process [10]. In fact, while 
some studies have supported the hypothesis [13-19], others have 
dismissed it [20-29]. Although most of the evidence relies on 
compositional pattern analysis at a regional scale [21], [30], 
recently some authors have implemented experimental approaches 
in communities of microorganisms, providing support to DNH 
[15], [17]. Even though the ecology of microorganisms is governed 
by processes equivalent to those that occur in multicellular 
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organisms [31-32], the specific mechanisms that promote 
successful invasion can differ considerably [32]. 

In the present study we evaluated DNH in experimental 

assemblages constituted by vascular plants. For this purpose, plants 
belonging to a wide taxonomic spectrum were used to estabhsh 
recipient experimental assemblages, which were then inoculated 
with seeds of a colonizing plant [Lactuca sativa, hereafter Lactuca). 
The experiments were composed of species that differed in their 
degree of phylogenetic relatedness with respect to Lactuca, a fact 
that allowed us to assess the effect of that factor on the colonizing 
success of the inoculated spci i(;s. As far as we can tell, this is the 
first time that DNH is evaluated experimentally in multicellular 
communities, specifically in plants. 

Methodology 

Experimental design 

Our experiment involved a total of 15 plant species (Table 1), 
14 of which were used to establish receiver assemblages and one 
[Lactuca) was used as colonizing or invading species of those 
assemblages. Five assemblages were formed and each of them was 
made up of a subset of five species of the 1 4 that were available 
(Table 1). Initially, these assemblages were organized in a 
taxonomic gradient from strong to weak relatedness with respect 
to Lactuca, and this was later confirmed by means of evolutionary 
distance metrics (see below). The experimental assemblages were 
designated Al, A2, A3, A4 and A5 (Table 1), and each of them 
was replicated eight times. For a control treatment (C), eight plots 
were set up to evaluate the colonization of Lactuca in a 
monoculture regimen. 

The assemblages were set up in small wooden plots 
(1x1x0.3 m) located in a greenhouse. A homogeneous mixture 
of sand and compost (1:1) was used as substrate, and evenly 
distributed in the plots. The main herbivores above ground (i.e., 
insects, birds and mammals) were completely excluded from the 
greenhouse. The plots were arranged in a pattern of equidistant 
rows and columns 2 m apart numbered sequentially. A random 
procedure allowed assigning each plot to a given kind of 
assemblage. In turn, each plot was subdivided into a grid of 
5x5 cm cells which were numbered consecutively from 1 to 400, 
and the species that would go in each cell in particular were 
assigned randomly. Seeds were planted in excess in each cell in 
order to ensure the future presence of plants, and once the 
seedlings had become established in each cell, the excess was 
removed by careful cutting with shears, so that aU the species had 
the same abundance in each plot (i.e., 66 cells per plot, occupying 
a total of 330 cells per plot), leaving 70 cells empty for the later 
planting of the colonizing species. 

The plots were watered periodically by means of a semi- 
automated system, in this way ensuring a homogeneous availabil- 
ity of water to the plants. This system consisted of an arrangement 
of rotary Micro-Jet sprinklers placed at a height of 1 .5 m; they 
were arranged in rows and columns equidistant from each plot to 
ensure uniform watering. Watering was performed at field 
capacity every three days. 

Three months after starting the planting in the assemblages (the 
seedlings had reached a height SlO cm), the Lactuca seeds were 
sown emulating the invasion of an already established receiver 
community. This planting process was carried out simultaneously 
in all the experimental assemblages and in the control treatment 
(from C, A1-A5) at a density of 70 seeds per plot. 



Lactuca colonization success 

The colonization success oi Lactuca was measured by means of 
five indicators. First, germination was evaluated by quantifying the 
number and percentage of cells with germinating seeds in each 
plot. Due to the fast germination of Lactuca, this indicator was 
measured three weeks after sowing the seeds (May 2010), which 
were considered to have germinated when their epicotyl had 
grown &2 cm. Second, the indicator of colonization success was 
the survival of Lactuca plants, and for that purpose, toward the 
end of the experiment [January 201 1) the number and percentage 
of surviving plants of the cohort existing at the beginning of the 
experiment were recorded. The third indicator considered the 
ratio of plants that reached the flowering stage, counting the 
number of plants that had one or more flowers in relation to the 
total number of Lactuca plants present in the plot. The fourth 
indicator evaluated the growth (height above the ground; cm) 
achieved by a random sample of 20 plants in each plot in January 
2011. The fifth indicator was recruiting, counting the number of 
new Lactuca seedlings that appeared spontaneously in the plots 
(number per plot). 

In order to control the concomitant effect of growth of the 
different species present on the studied assemblages, the height 
achieved by the plants in each assemblage in each plot was 
measured. For that purpose a random sample of 20 plants was 
taken, measuring twice their height (cm) above the ground: at the 
beginning of the experiment (May 2010) to evaluate their effect on 
the germination ot Lactuca seeds, and at the end of the experiment 
(January 2011) to evaluate their effect on the growth of Lactuca 
plants. 

Phylogenetic relatedness among species 

A metaphylogeny was reconstructed for the 15 species of 
angiosperms included in our study using Phylomatic [33] version 
R20031202, which is based on the APG III phylogeny [34]. 
Within the family Apiaceae, the phylogenetic relations between 
Anethum, Coriandrun and Petroselinun were resolved using [35] . 
The topology of the resulting tree was cahbrated by age, based on 
the divergence times documented by [36], using the BLADJ 
algorithm implemented in Phylocom [37]. 

Based on this calibrated tree, two phylogenetic diversity indices 
were calculated to characterize- (-acli assemblage. On the one 
hand, the average length of the branches was calculated for aU the 
pairs of species of the assemblage (mean phylogenetic distance, 
MPD; [38-40]) both before (MPD^„) and after (MPD^,„,;) planting 
Lactuca. Furthermore, the average length of the branches between 
every species and its nearest neighbor in the assemblages was 
CEilculated (mean nearest neighbor distance, MNND; [38], [40]) 
both before (MNND^^) and after (MNND^oj^) planting Lactuca. 
The values of MPD and MNND are expressed in millions of years 
(my). Because there was no loss of species over the course of the 
experiments, the post-invasion phylogenetic diversity shows the 
effect of adding Lactuca to the plots. 

As indicators of phylogenetic relatedness between Lactuca and 
the remaining species in each assemblage, we used the average of 
MPD between Lactuca and each member of the assemblages 
(namc-d MPD/^„ri„„,_), and the distance bet\\x'cn Lactuca and the 
nearest neighbor within assemblages (named MNND/,„,.(,„,.„). The 
two metrics provide different ways in which phyiogenetic 
relatedness can be conceived in the particular DNH context, 
because ^iPG Lactuca Considers values at the community level, while 
MNND^actoca includes only the nearest neighbor. 
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Table 1. Composition of the experimental assemblages and characterization of their phylogenetic diversity based on evolutionary 
distances. 



Assemblages (Treatments) 



Plant species 


A1 


A2 


A3 


A4 


A5 


Asteraceae 


Lactuca sativa L. 


1 


1 


1 


1 


1 


Matricaria ciiamomilla L. 


1 


0 


0 


0 


0 


Apiaceae 


Anettium graveolens L. 


1 


1 


0 


0 


0 


Petrosetinum crispurri (Mill.) Fuss 


1 


1 


0 


0 


0 


Coriandrum sativum L. 


1 


1 


0 


0 


0 


Solanaceae 


Capsicum baccatum Jacq. 


0 


1 


1 


0 


0 


Solarium metongena L . 


0 


1 


1 


0 


0 


Lamiaceae 


Ocinum basilicum L. 


0 


0 


1 


0 


0 


Brassicaceae 


Brassica oleraceae L. 


0 


0 


1 




0 


Eruca sativa Mill. 


0 


0 


1 




0 


Fabaceae 


Pisum sativum L. 


0 


0 


0 






Trifotium repens L. 


1 


0 


0 






Vicia atropurpurea Desf. 


0 


0 


0 






Vicia faba L. 


0 


0 


0 






Poaceae 


Zea mays L. 


0 


0 


0 






MPDp„ (my) 


182.3 


165.1 


208.2 


142.6 


184.8 


MPDp„„ (my) 


187.1 


184.0 


219.4 


179.7 


212.2 


MNNDp„ (my) 


120.0 


59.2 


94.8 


64.0 


109.2 


MNNDp„„ (my) 


141.8 


142.6 


143.8 


143.0 


143.3 


MPDi^mcc (my) 


196.8 


222.0 


242.0 


254.0 


267.6 


MNNDt„„„„ (my) 


88.0 


214.0 


234.0 


254.0 


254.0 



MPDp^e and MPDposf represent the average distance of the branch lengths between the pairs of species of the assemblage before and after the invasion by Lactuca, 
respectively; MNNDp^p and MNNDpo^j represent the average distance of the branch lengths between each species and its nearest neighbor in the assemblage; MPDtoctuco 
corresponds to the average branch length between Lactuca and each of the members of the assemblage, respectively; and MNNDtacfuca was calculated as the branch 
length between Lactuca and its nearest neighbor, my: million of years. 
doi:1 0.1 371 /journal.pone.Ol 05535.t001 



Statistical analyses 

In a first group of analyses we compared the colonization 
success indicators ot Lactuca in monocultures (control treatment, 
C) versus those recorded in the assemblages (treatments A1-A5), 
using one-way analysis of variance (ANOVA) for germination, 
survival, flowering, growth, and colonization of Lactuca in the 
different treatments. Prior to each analysis, all the dependent 
variables were transformed into logarithmic functions (In (x)). The 
Tukey test was apphed to recognize those treatments that showed 
statistical significance. 

With the purpose of assessing the relation between phylogenetic 
distance and the colonization success of Lactuca, we performed 
covariance analyses (ANCOVA) for the different assemblages. 
Here the MPD/,„c(„„, and MNND j^tiua distances were concom- 
itant factors, while the colonization indicators were considered 
dependent variables. In these analyses, the average as well as the 
variance of the height reached by the resident plants in the 



assemblages were included as covariables. AH these factors were 
normalized by means of the logarithmic function (In (x)). In all our 
analyses we used the Type III sum of squares. 

Results 

Pre- and post-invasion phylogenetic diversity 

Before the inoculation with Lactuca, the assemblages showed 
MFDpre values varying between 142.6 and 208.2 my, and 
MNND^„ values between 59.2 and 120.0 my (Table 1). After 
the inoculation with Lactuca, the MPDp„^t values increased to 
between 179.7 and 212.2 my, while those of MNNDyj„,; were 
between 141.8 and 143.8 my (Table 1). Therefore, the inoculation 
of Lactuca increased significantly the phylogenetic diversity values 
measured as MPD and MNND (Table 1; Wilcoxon Signed Ranks 
test; in both cases T= 7.5; P<0.05). 
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Colonization success 

The average germination rate of Lactuca among the treatments 
varied between 55.8 and 98.3% (Figure lA), with statistical 
differences between them {F-bA.2; d.f. = 5; P<0.05). These 
differences were determined by the greater germination in the 
control treatment with respect to the experimental assemblages, 
which did not show differences between them (Figure lA). 

The average growth of Lactuca in the different treatments 
varied between 24 and 60 cm (Figure IB). These values showed 
statistical differences among treatments (F — 5.8; d.f. — 5; P<0.05), 
determined by the control treatment and the A2 assemblage with 
respect to the other assemblages (Figure IB). On the other hand, 
the average survival ot Lactuca varied between 56.8 and 93.8% 
(Figur(; IC), \-alucs that also showed significant differences among 
compared treatments {F=5.7; d.f. —5; P<0.05); again these 
differences were determined by the control treatment with respect 
to the other assemblages (Figure IC). 

With respect to flowering, this measurement varied between 9.0 
and 33.6% (Figure ID), with statistical differences among the 
treatments {F= 10.4; d.f. = 5; P<0.05), which were established by 
the control treatment with respect to the remaining species in the 
assemblages (Figure ID). Finally, the average recruitment varied 
between 1.3 and 8.5 seedlings per plot (Figure IE), showing 
statistical differences among the treatments {F = 10.4; d.f. = 5, P< 
0.05), which were determined by a higher seedling density in the 
control treatment with respect to that of the experimental 
assemblages (Figure IE). 

DNH 

The phylogenetic relatedness between Lactuca and the resident 
members in eac:h of the assemblages varied between 196.8 and 
267.6 my for MPD^^^ju^^ and between 88.0 and 254 my for 
M!<n<IDLa,tuca (Table 1). The ANCOVA results showed that 
neither MPD/ nor MNNDiari„„, had significant effect on 
the germination, growth, survival, flowering, or recruitment of 
Lactuca (see Table 2), showing that the fi\c indicators ot Lactuca 
colonization success were not affected by the phylogenetic 
relatedness between this species and the experimental assemblages. 
The MNND^o,.j„,,j distance showed a marginally significant effect 
on the growth and recruitment of Lacluca (p = 0,052 and 0.085, 
respectively; see Table 2). However, these results were due to 
effects recorded in A2 and A4 treatments, no clear trend in 
support (or refutation) of DNH. 

The average height of the assemblages had a significant effect 
on three of the colonization success indicators, specifically on 
germination, growth, and flowering of Lactuca (Table 2). This 
effect was determined by a decrease in the colonization success on 
plots with taUer plants. Finally, the variance of the height of the 
assemblages did not have significant effects on any indicator of the 
establishment of Lactuca (Table 2). 

In brief, the colonization success of Laclnca, measured by 
means of five different indicators (germination, growth, survival, 
flowering, and recruitment) was significandy greater in the control 
treatments (i.e., Lactuca monocultures) than in the experimental 
assemblages, showing that in the presence of other species the 
colonization success of Lactuca is reduced. However, Laclnca 
colonization did not show significant differences between the 
experimental assemblages, regardless of the phylogenetic related- 
ness with the receiver assemblage (measured as MPD/^rf^^o and 
MNND/ „,.(,„„). Finally, the average height achieved by the plants 
of each assemblage reduced three of the five indicators of 
colonization success by Lactuca. 



Discussion 

DNH states that if the competitive interactions between 
phylogeneticaUy close species are more intense, so the colonization 
success wiU be reduced when a given taxon colonizes communities 

that contain related species. Conversely, the colonization success 
will increase if the invasion occurs in communities consisting of 
phylogeneticaUy distant taxa. Using an experimental gradient of 
phylogenetic relatedness between five plant species (receiver 
assemblages) and a colonizing species (Lactuca sativa), we found 
that phylogenetic relatedness did not influence the colonization 
success of the inoculated species. Therefore, these findings do not 
support DNH, and this is valid for a combination of five measures 
of colonizing success (germination, growth, survival, flowering, and 
recruitment) and two of phylogenetic relatedness (MFDj^tMa and 
MNNDi,,,,„,„). 

From the perspecti\'e of DNH, our experiment shows two 
important limitations. First, the experimental plots did not cover 
the (continuous) spectrum of phylogenetic relatedness between 
colonizer and recipient assemblage. However, this does not seem 
to have affected our analyses because using a posteriori tests we 
found no differences between Lactuca colonization success 
recorded in treatment Al (MPDi„rf„,„ = 196 my; MNNDi,,d„,„ = 
88 my) and the other treatments (MPT)i^^i„ca — 222-267 my; 
MNNDiarf„,o = 214— 254 my). Second, our experimental design 
implied that not only the phylogenetic relatedness means greater 
competitive intensity, but this can also be due to the transfer of 
specialist natural enemies (e.g., herbivores and parasites) from the 
receiver assemblage to colonizing species [41,42]. This also 
reduces the establishment success of colonizing species, a factor 
that was not assessed in our experiments because we exclude the 
effect of herbivore interactions. Moreover, the five indicators of 
successful colonization ot Lactuca showed inhibitory effects in the 
experimental assemblages compared to the control treatments 
[Lactuca monocultures). In mechanistic terms, these results 
probably reflect inhibitory interactions that restrict access to light 
and nutrients, exercised by species from recipient assemblages 
upon Lactuca. 

Our results differ with respect to two studies that so far have 
evaluated DNH experimentally. Among previous experimental 
work, one study [15] used a Serratia marcescens strain as 
colonizing taxon, while another strain of the same species was 
part of the receiver communities; and another, [17] used 
colonizing species of Metschnikowia and Candida, which had 
congeneric representatives in the recipient assemblages. In 
contrast to these studies, in our experiments we worked with 
more distant levels of relatedness, since between Lactuca and 
Matricaria chamomilla, both species belonging to the same 
taxonomic family (Asteraceae), we found an MNND distance of 
88.0 my. Therefore, the use of the congeneric or same strain 
(conspecific) taxa may magnify the effect of competition between 
colonizer and receiver assemblages [41], explaining —at least 
partially- our discrepant results. 

A hypothesis opposed to DNH is related with preadaptation 
[29], [42-43], which proposes that phylogenetic closeness 
promotes the colonization process through facilitation among 
related taxa [41-42], [44-45]. Our results do not support this 
hypothesis, because colonization by Lactuca did not show a 
positive association with the phylogenetic relatedness of the 
assemblages. Moreover, the most important mechanisms of 
Lactuca colonization were related to the average height of the 
plants of the receiver assemblages. This trait, which is an indicator 
of the availability of light to the colonizer affected the germination, 
growth, and flowering response of Lactuca. 
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A B 




C A1 A2 A3 A4 A5 C A1 A2 A3 A4 A5 



C D 




C A1 A2 A3 A4 A5 C A1 A2 A3 A4 A5 



E 




C A1 A2 A3 A4 A5 



Figure 1. Colonization success (average ± S.D.) of Lactuca invading experimental plant communities. Colonization success was assessed 
as germination (%) (panel A), growth (cm) (panel B), survival (%) (panel C), flowering (%) (panel D), and recruitment (seedlings per m'^) (panel E) of 
inoculation of Lactuca recorded in monoculture (control, C) and five assemblages {A1-A5) with which it has different degrees of phylogenetic 

relatedness. 

doi:1 0.1 371/journal.pone.01 05535.g001 
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Our results suggest that under the spatial scale used there was a 
nuU effect of the phylogenetic distance on colonization success and 
the invasion's contingent conditions such as the current properties 
of resident assemblage, and the intrinsic properties of invaders 
would be more relevant in determining colonization success than 
factors associated with phylogenetic relatedness, as has been 
discussed by other authors [10], [41], [43], [46], [47]. This does 
not mean that phylogenetic distance between invader and resident 
assemblage does not affect the colonization process. Indeed, the 
phylogenetic distance would have a non-linear effect on coloni- 
zation process [43], which would require increasing the spatial 
scale of the studies, increasing the number of species in resident 
assemblages, or performing comparative studies that consider a 
larger number of invader species in several recipient communities 
[41]. In the absence of more evidence and experimental contrasts, 
it seems premature to accept the null effect as a definitive answer. 
However, the available evidence at least allows questioning the 
generality of DNH, since it has been verified experimentally only 
when the colonizing species are closely related to the members of 
the receiving community [15,17], while if the colonizer is 
phylogeneticaUy more distant (this study), the colonization success 
becomes independent of the evolutionary relatedness. A recent 
study [21] evaluated the change in the intensity of the interactions 
between vascular plant species along a phylogenetic gradient (with 
MNND distances ranging between 0 my and 81 my). Although 
these experiments were made between pairs of species and not in 
communities, the support of DNH was only partial because the 
evolutionary closeness not only allowed greater inhibitory inter- 
action intensity (such as competition), but it also increased the 
intensity of facilitating interactions, which have an effect that is the 
opposite of that expected by DNH [10,43]. From the community 
viewpoint it would be extremely important to quantify the 
combined effect of the inhibitory and facUitative interactions on 
the success of colonization, and how this balance is expressed 
along the colonizer-community relatedness gradient. Along this 
line, the evolutionary distance metrics should include most of the 
components of the community (i.e., herbi\ ores, parasite's, facilita- 
tors) with which a colonizer can interact, and not only the 
taxonomical composition of assemblages. These kinds of efforts 
can be difficult to implement under field conditions, but 
experimental or modeling approximations can assist in disentan- 
gling this complexity [15]. We agree with two previous studies 
[10,43], which proposed deeper studies that allow establishing the 
role of the phylogenetic structure of the communities in their 
susceptibility to being invaded. 



An important challenge that still needs to be elucidated is how 
these experimental conditions reflect the heterogeneity and 
complexity of the invasive processes under field conditions. For 

example, it is feasible that invasion by congeneric or conspecific 
species occurs when a taxon expands its range by means of a 
reaction-diffusion process [48], colonizing adjacent communities, 
which likely contain closely related taxa. In contrast, invasion that 
involves spread to new continents or distant regions probably 
represents better the case in which the colonizing species has a 
rather distant relationship with the members of the receiving 
community [48]. 

In summary, our results do not show an association (positive or 
negative) between phylogenetic distance and the colonizing success 
of an inoculated species, so they do not support DNH. In view of 
the small number of studies that have evaluated DNH, particularly 
from the experimental standpoint, it may be premature to 
generalize on the role of phylogenetic relatedness in determining 
the result of the invasive process. However, our results and the 
a\ ailable background information suggest that DNH may have an 
explanatory domain restricted to only one part of the phylogenetic 
spectrum. 
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Table SI Colonization success of Lactuca recorded in 
each experimental plot. The colonization indicators were: 
Germination (%), Growth (cm). Survival (%), Flowering (%) and 
Recruitment (No./m2). A1-A5 represent different treatments of 
recipient plant assemblages and C the control treatment (i.e., 
Lactuca in monoculture). MPD/^ctea MNND/^ftjca are two 
metrics of phylogenetic distances for Lactuca and the recipient 
assemblages. MPDi^jy^a is the average distance between Lactuca 
and each member of the assemblages; MNND£,jrf„,.„ is the distance 
between Lactuca and its nearest neighbor in each assemblage. As a 
concomitant factor, the Height (average and variance; both in cm) 
achieved by the plants in each plot at the end of the experiment 
was considered. 
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